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Abstract 

Using the integral equations of the Noncrossing Approximation, the differ- 
ential conductance is computed as a function of voltage for scattering from a 
two channel Kondo impurity in a point contact. The results compare well to 
experimental data on Cu point contacts by Ralph and Buhrman. They sup- 
port a recently proposed scaling hypothesis, and also show finite temperature 
corrections to scaling in agreement with experiment. The conductance signal 
is predicted to be asymmetric in the bias when the impurity is not equally 

coupled to left and right moving electrons. 

PACS numbers: 72.10-d, 72.15.Fk, 72.10.Qm, 63.50.+X 
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The two channel Kondo model |T|, or equivalently the Kondo limit of the two channel 
Anderson model has been applied to a wide variety of interesting physical systems: 
heavy fermion compounds high T c superconductors |J, and two level systems (TLS) in 
metals 0-0- Although this model contains many of the salient features of the experiments, 
e.g. marginal- Fermi-liquid behavior [ID|,|5|, the experimental proof for the physical existence 



of two channel Kondo impurities is far from certain. 

One of the strongest experimental cases for the existence of two channel Kondo behavior 
is an experiment by Ralph and Buhrman |§ on clean Cu point contacts. At low temperatures 
their samples show the correct temperature, magnetic field, and voltage dependence for a two 
channel Kondo impurity such as a TLS with electron assisted tunneling. Recently ||, the 
data has also been shown to be consistent with a scaling ansatz motivated by the equilibrium 



Conformal Field Theory (CFT) solution of the problem [llj]. In this paper we compute the 
differential conductance for a two channel Kondo impurity in a point contact and find 
excellent quantitative agreement with the above experiment, lending strong evidence for the 
existence of two channel Kondo impurities in this system. 

The two channel Kondo model consists of two kinds or channels of mutually noninter- 
acting conduction electrons which are coupled to the same impurity spin via an exchange 
interaction jj]]. For the physical realization of interest here |§ a TLS plays the role of the 
impurity spin. The states of the TLS correspond to the impurity spin-up and spin-down 
states. The electrons scattering from the impurity may be characterized by parity. The 
different parity states play the role of spin-up and spin-down electrons in the conventional 
Kondo problem. Thus, parity plays the role of the active degree of freedom altered by the 
interaction, while the physical spin is a spectator degree of freedom, which is conserved by 
the interaction and which constitutes the two channels of the model |J. 

It is convenient to represent this system by an Anderson hamiltonian of a particle on an 
impurity level far below the Fermi surface hybridizing with the two channels of conduction 
electrons. We compute the differential conductance within the Non Crossing Approximation 



(NCA) for the infinite U Anderson model in the Kondo limit ||T2Hl^1. The NCA has been 



very successful in describing the one channel Kondo problem except for the appearance 
of spurious nonanalytic behavior at a temperature far below the Kondo temperature T#. 
These spurious low-T properties are due to the fact that the NCA neglects vertex corrections 
responsible for restoring the low temperature Fermi liquid behavior of the one channel model 
[Tj|. However, it has recently been shown [T5J that for the two channel problem, where the 



complications of the appearance of a Fermi liquid fixed point are not present, the NCA does 
give the exact low-frequency power law behavior of the impurity spectral function Ad(u)) 
at zero temperature. Therefore, we expect to achieve a correct description for quantities 
involving A d (like the conductance). 

In the NCA approach the electron on the impurity is represented in terms of a pseudo 
Fermion operator, /, and a slave Boson operator, b fl7| . In order to calculate the conductance 



at finite bias, the NCA must be generalized using nonequilibrium Green functions [18]. One 



solves for both the retarded Green functions for these operators, G{ and G b , and for the 
'lesser' Green functions G< and G b K , which contain information about the nonequilibrium 
distribution function. The integral equations for these four functions have been solved for the 
one channel case by Meir, Wingreen, and Lee [ITUJ. However, our numerical implementation 



pi| is different from theirs, and we are able to go more than two orders of magnitude lower 
in temperature, deep into the low temperature scaling regime described below. 

In the Anderson model the characteristic energy scale, the Kondo temperature, Tk, 



depends on the width, T, and position, E d <C — T, of the bare Anderson impurity level |L3 
Besides the voltage, V, and temperature, T, we can also vary the couplings of the impurity 
to left and right moving electrons and T^, normalized to the width T (T^ + = 1). 
Except for Fig. 4 all the numerical results have Tl = T^. 

Letting F eff (u) = T L F(u+eV/2)+T R F(u-eV/2), where F(u) = l/(l + e^), and using 
the conventions of Miiller-Hartmann |TJ] for the spectral functions, A(u) = —ImGl{uj) /it, 
B(uj) = — ImGr(u))/7r, and the 'lesser' Green functions, a (a?) = G<(u;)/27r and b(ui) = 
G b < (u)/2ir, the nonequilbrium NCA equations for the N = 2 spin, M = 2 channel Anderson 
model are [^l[ 



Jg^nv/^ + ^w (i) 
KlSr rM /7 B <— ni-F, ff i,)} (2) 
V%fa = ™jT« u + W-W<)\ (3) 

The true impurity spectral function, A^, is computed from the slave Green functions via the 
convolution 

A d (u) = [-[ a(e)B(e - u) + A(e)b(e - u)) . (5) 

J 7T 

A point contact consists of two leads joined by a small constriction. Any additional 
scattering in the vicinity of the constriction should cause a decrease in the conductance 
because it impedes the flow of electrons. On the other hand, in a tunnel junction, tun- 
neling may be assisted by electrons hopping on impurities in the junction, increasing the 
conductance. Thus, it is not surprising that when we generalize earlier calculations for the 



nonlinear current through a tunnel junction f22|,|23[ to the case of a point contact, we find a 



similar expression for the current, except for an overall minus sign: 

I- J oc - J duA d (u)[F(u - V/2) - F(oj + V/2)], (6) 

where I and Iq are the currents with and without the impurity. In deriving Eq. ^| we have 
assumed that the point contact is clean, namely the transmission coefficients are close to 
unity, and that all hopping matrix elements are slowly varying on the scale of the Kondo 
temperature. The conductance, G(V,T), is computed by taking a numerical derivative, 
dl/dV. We will assume that the background conductance dlo/dV is ohmic. 

In Figure 1 we show the zero bias conductance, G(0, T), computed in this manner (T^ = 
T R ). As expected [II]], the conductance shows a T 1//2 dependence at low temperatures with 
deviations starting at about 1/4 T%. The Kondo temperature is determined by the width 
at half maximum of the zero bias impurity spectral function, Ad, at the lowest calculated 
temperatures (see inset). The slope of the T 1//2 behavior defines the constant B^,: 



G(0,T)-G(0,0)=B E T 1 / 2 . (7) 

The experimental data also show a T 1 / 2 - dependence, but it is difficult to deduce an accurate 
estimate of T% by looking at the deviations from T 1 / 2 behavior (assuming that they would 
occur at 1/4 Tk). An educated guess gives Tk to be about 8 Kelvin for sample 1 and 2 and 
significantly less for sample 3 of Ref. |J. 

Recently, it has been proposed from a CFT solution of the problem in equilibrium that 
the experimental data show scaling of the conductance G as a function of voltage bias V 
and temperature T of the form || 

pV 

G(V, T) - G(0, T) = BxT^HdA—)), (8) 

Kb± 

where if is a universal scaling function (H(0) = and H(x) ~ x 1 ^ 2 for x >>1) and B^ and 
A are nonuniversal constants. In order to examine whether this ansatz is correct, in Fig. 
2 the rescaled conductance is plotted as a function of (ey/Zc^T) 1 / 2 for the numerical data 
(a) and the experimental data (b, for the best sample (1)). Motivated by the symmetry 
of the experimental conductance-voltage curves we choose equal coupling to left and right 
moving electrons, I\ = (see also Fig. 4). Considering that after fixing B-% using Eq. 
(7) there are no adjustable parameters, the agreement is extraordinary. The collapse of the 
various T curves at low bias and the linear behavior for the low T curves in the range of 
2 < (eV/fceT) 1 / 2 < 4 is in agreement with the proposed scaling ansatz Eq. |8|. However, the 
slope of the linear part shows temperature dependence for both the experimental and the 
numerical data. This is not contradictory to the scaling ansatz, but it does show that there 
are significant temperature dependent corrections to scaling within the present approach. 

To analyze the scaling plots in more detail, the low bias portion is replotted in Fig. 
3(a). The conductance follows an approximate (eV/ksT) 2 - behavior even for (eV/ksT) > 1 
before it levels off and enters the (eV/ksT) 1 ^ 2 - region at higher bias. The prefactor of the 
quadratic dependence shows no observable temperature dependence until approximately 
0.1T K and consequently obeys the scaling ansatz. 
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In Fig. 3(b) the slopes of the straight line fits of the linear regions in Fig. 2 are 
plotted as a function of temperature. Both the numerical data and the experiment show 
clear temperature dependence. Although By, may be determined directly from the zero 
bias conductance, the Kondo temperature is more difficult to determine experimentally. In 
Fig. 3(b) we have chosen values for the Kondo temperature which are consistent with the 
estimates from the deviation of the zero bias conductance from T 1//2 behavior. The resulting 
curves are in good quantitative agreement. In order to show that the experimental and 
numerical curves indeed coincide (for a given T/Tk) we also compare the intercepts of the 
straight line fits for the same Kondo temperatures (inset Fig. 3 (b)). The numerical data 
do fall right in the middle of the scatter from the three different samples ||24|| . Note that 



our theory does not have the additional parameter, A, of Eq. ||, but adjusting the Kondo 
temperature has some of the same effect as adjusting A. 

Finally, we have also studied the conductance in cases where the impurity is asymmetri- 
cally coupled to left and right moving electrons. The NCA equations have no symmetry for 

<-> T/j; however, they are symmetric under the parity operation «-> T r and V «-> —V. 
Figure 4 shows the conductance signal one should observe in an asymmetrically coupled 
sample. The conductance is clearly a nonuniversal function of the couplings and Yr. 
The experimental data are symmetric to within experimental error, indicating that left and 
right moving electrons are equally coupled to the TLS. We think this is reasonable for a 
point contact where the impurity has no particular asymmetry with respect to left and 
right moving electrons, but not for a tunnel junction, where one would expect exponential 
dependence of the couplings on the location of the impurity. 

In conclusion, we have performed numerical evaluations of the NCA integral equations 
for the two channel Anderson model out of equilibrium. We find outstanding agreement 
with data from an experiment on Cu point contacts ||. Scaling of the conductance at low 
bias (eV < fc^T) and temperatures is verified. As V and T are increased, the calculation 
exhibits finite-T corrections to scaling, again in agreement with experiment. Asymmetric 
conductance signals are predicted for impurities asymmetrically coupled to the leads. Thus, 
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this work lends strong support to the existence of two channel Kondo impurities in the Cu 
point contacts of Ralph and Buhrman. 
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(M.H. and S.H.), and by a Feodor Lynen Fellowship of the Alexander v. Humboldt Found. 
(J.K). 



7 



REFERENCES 

[1] P. Nozieres and A. Blandin, J. Phys. (Paris) 41, 193 (1980). 
[2] J. R. Schrieffer and P. A. Wolff, Phys. Rev. B149, 491 (1966). 

[3] D.L. Cox, Phys. Rev. Lett. 59, 1240 (1987); Physica (Amsterdam) 153-155C, 1642 
(1987); J. Magn. Magn. 76&77, 53 (1988). 

[4] C. L. Seaman et. al, Phys. Rev. Lett. 67, 2882 (1991), B. Andraka and A. M. Tsvelik, 
Phys. Rev. Lett. 67 , 2886 (1991). 

[5] D. L. Cox et. al, Phys. Rev. Lett. 62, 2188 (1989), V. J. Emery and S. Kivelson, 
Phys. Rev. B 46, 10812 (1992); Phys. Rev. Lett. 71, 3701 (1993), T. Giamarchi et al, 
Phys. Rev. Lett. 70, 3967 (1993). 

[6] A. Zawadowski, Phys. Rev. Lett. 45, 211 (1980; K. Vladar and A. Zawadowski, 
Phys. Rev. B 28, 1564, 1582, 1596 (1983). 

[7] A. Muramatsu and F. Guinea, Phys. Rev. Lett. 57, 2337 (1986), S. Katayama et. al, 
J. Phys. Soc. Jpn. 56, 697 (1987), T. Ishiguro et. al, Phys. Rev. Lett. 69, 660 (1992). 

[8] D. C. Ralph and R. A. Buhrman, Phys. Rev. Lett. 69, 2118 (1992). 

[9] D. C. Ralph, et. al. , Phys. Rev. Lett. 72, 1064 (1994). 

[10] C. M. Varma, et. al. , Phys. Rev. Lett. 63, 1996 (1989). 

[11] A. A. W. Ludwig and I. Affleck, Phys. Rev. Lett. 67, 3160 (1991); I. Affleck and 
A. A. W. Ludwig, Phys. Rev. B 48, 7297 (1993). 

[12] P. Coleman, Phys. Rev. B 29, 3035 (1984). 

[13] E. Muller-Hartmann, Z. Phys. B 57, 281 (1984). 

[14] N. E. Bickers, Rev. Mod. Phys. 59, 845 (1987). 

[15] T. A. Costi et. al. , submitted to Phys. Rev. Lett. (1994). 

8 



[16] D. L. Cox and A. E. Ruckenstein, Phys. Rev. Lett. 71, 1613 (1993). 
[17] S. E. Barnes, J. Phys. F6, 1375 (1976); F7, 2637 (1977). 

[18] D. C. Langreth, 1975 NATO Advanced Study Institute on Linear and Nonlinear Electron 
Transport in Solids, Antwerpen, 1975 (Plenum, New York, 1976), B17, p. 3. . 

[19] Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 70, 2601 (1993). 

[20] For the equilibrium case see J. Kroha, Ph.D. Thesis, U. Karlsruhe, (1993). 

[21] Instead of the flat density of states (DOS) of Eqs. (1-4), we use a gaussian DOS, which 
can be included in F e ff. 

[22] S. Hershfield, J. H. Davies and J. W. Wilkins, Phys. Rev. Lett. 67, 3270 (1991). 

[23] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 (1992). 

[24] We can also directly identify the experimental curves with the numerical data over the 
scaling range. 



9 



FIGURES 

FIG. 1. Temperature dependence of the zero bias conductance (Tl = Tr). The zero bias 
conductance has T 1 / 2 - dependence for T < T&-/4. This can be used to roughly extract Tk from 
the experimental data. Inset: The impurity spectral function Ad(oj) for several voltages. The 
width at half maximum of the zero bias spectral function (dotted curve) determines Tk- As the 
voltage is increased to eV = ksTx the Kondo resonance is reduced (solid curve). At very large 
bias the resonance shows a shoulder and eventually two peaks. In this paper we compare theory 
and experiment in the scaling regime, T <C Tk- 

FIG. 2. Scaling plots of the conductance for (a) theory and (b) experiment [13]. With Tl = Tr 
and .Be determined from the zero bias conductance (see Fig. 1), there are no adjustable parameters. 
There are roughly two regimes in these plots. For (eV/ZcsT) 1 / 2 < 1.5 the curves collapse onto a 
single curve and the rescaled conductance is proportional to (eV/kBT) 2 . For 2 < (eV/ksT) 1 / 2 < 4 
the rescaled conductance is linear on these plots. There are substantial corrections to scaling even 
at temperatures small compared to Tk (see figure 3 (b)). At even larger biases this linear behavior 
rounds off, indicating the breakdown of scaling. The temperatures in the theory and experimental 
plots are in units of Tk and Kelvin, respectively. 

FIG. 3. Quantitative analysis of the scaling plots at (a) low bias and (b) high bias, (a) The 
low bias rescaled conductance as a function of (eV/ksT) 2 (curves are offset). Both theory and 
experiment (sample 1 of Ref. [9]) show quadratic behavior at low bias. The symbols correspond 
to the temperatures shown in Fig. 2 . (b) The slope of the straight line fits of the linear region 
in Fig. 2 as a function of temperature for sample 1 (o,+), sample 2 (<>,□), sample 3 (x,*) (for 
V > 0, V < 0) and the numerical data (A). The Kondo temperature for the experimental curves 
is consistent with the deviation of the zero bias conductance from T 1 / 2 - behavior. All curves drop 
with increasing temperature even at temperatures small compared to the Kondo temperature. The 
inset shows the behavior of the intercepts of the fits. 
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FIG. 4. Dependence of the conductance on the coupling of the impurity to left (right) moving 
electrons, T L m\, with Tl + Tr = 1. In Figs. 1-3 Tl equals Tr. For a typical low temperature 
(T = O.OlTft;) one observes an asymmetric signal, G(V, T) / G(— V, T) for Tl ^ Tr. Note that B^,, 
which depends on Tl, has been divided out. These curves show that the scaled conductance is not 
universal if one allows for asymmetric coupling of the impurity to left and right moving electrons. 
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